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Abstract. An overview is given of multiwavelength observations of young supernova remnants, with a focus on the 
observational signatures of efficient cosmic ray acceleration. Some of the effects that may be attributed to efficient cosmic 
ray acceleration are the radial magnetic fields in young supernova remnants, magnetic field amplification as determined with 
X-ray imaging spectroscopy, evidence for large post-shock compression factors, and low plasma temperatures, as measured 
with high resolution optical/UV/X-ray spectroscopy. Special emphasis is given to spectroscopy of post-shock plasma's, which 
offers an opportunity to directly measure the post-shock temperature. In the presence of efficient cosmic ray acceleration the 
post-shock temperatures are expected to be lower than according to standard equations for a strong shock. For a number of 
supernova remnants this seems indeed to be the case. 
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INTRODUCTION 

Young supernova remnants are interesting from many 
different points of view: they are important for study- 
ing the freshly synthesized supernova material, thereby 
giving clues about the supernova explosion itself |[T] |2l 
for some recent examples], and they are of interest for 
many physical processes associated with low density 
plasmas, for example coUisionless shock physics and 
non-equilibrium ionization. However, in keeping with 
the theme of this symposium, I will focus here on those 
multiwavelength aspects that reveal something about the 
cosmic ray acceleration properties of young supernova 
remnants (SNRs). 

For a long time SNRs have been considered the prime 
candidates for providing the bulk of the cosmic rays ob- 
served on earth in the energy range up to £ 3 x lO'^ eV 
(the "knee" in the cosmic ray spectrum), but possibly 
even up to lO'^ eV (the "ankle"). The supernova ener- 
getics and rates are sufficient to explain the energy den- 
sity of Galactic cosmic rays, but the problem was that 
the conditions in SNRs seemed not right for accelerating 
particle up to, or beyond, energies of lO'^ eV. In partic- 
ular, the magnetic field strengths and turbulence needed 
for fast acceleration, were considered unlikely to exist in 
SNRs |3|. 

Since the 1950-ies it is known that at least some par- 
ticles get accelerated up to relativistic energies, since 
SNRs turned out to be sources of synchrotron radio emis- 
sion, which indicated the presence of relativistic elec- 
trons [4 for an early discussion]. However, the electrons 
producing radio synchrotron emission have energies in 
the GeV range, far short of the 10 '^eV needed to explain 



the cosmic ray spectrum. 

From the 1950-ies to the mid 1990-ies not much ob- 
servational progress was made concerning cosmic ray ac- 
celeration by SNRs. Disappointingly, y-ray experiments 
like COS-B, or EGRET revealed little evidence for the 
presence of high energy cosmic rays inside SNRs, al- 
though EGRET did find a possible connection between 
7-ray sources and old SNRs near molecular clouds [SJ. 

Over the last decade the situation has changed con- 
siderably. X-ray observations have shown that a) elec- 
trons are accelerated to > 10 TeV energies in some 
SNRs 0, b) magnetic fields in SNRs are much larger 
(> 100 /jG) than the compressed average field in the in- 
terstellar medium (-- 5 flG) [e.g. 7 , 8 , 9 , 10 1. 

The presence of > 10 TeV electrons reveals that SNR 
shocks are capable of acceleration particles to very high 
energies, although it does not yet provide direct evidence 
that ions, which constitute the bulk of the cosmic rays 
observed on earth, are accelerated beyond 10^^ eV. 

However, the relatively high magnetic field at least 
indicates the conditions for accelerating particles up to 
10^^ eV are present. Moreover, the high magnetic fields 
probably arise through plasma wave generation by cos- 
mic rays themselves [e.g. TT]. High magnetic fields in 
SNRs, therefore, provide indirect evidence for efficient 
cosmic ray acceleration. 

Direct evidence that, indeed, ions are accelerated to 
high energies by young SNRs comes from TeV 7-ray 
observations. Several SNRs have now been detected by 
Cherenkov 7-ray telescopes |il2] [13] HH [15], but the 
interpretation of the results are hotly debated lT6l : is the 
7-ray emission caused by inverse Compton scattering of 
> 10 TeV electrons, or is it due to the decay of neutral 



pions created by energetic ion-ion collisions? 

Although this is an ongoing debate, we will likely 
make major advances in our understanding in the next 
decade. One of the reasons is the rapid developments 
in TeV astronomy, but multiwavelength data are neces- 
sary to interpret the TeV data (e.g. to estimate densities, 
magnetic fields), and as I will show multiwavelength ob- 
servation will help to reveal whether ions are efficiently 
accelerated, where and when they will are accelerated, 
and how much of the shock energy will be taken up by 
cosmic rays. 

SUPERNOVA REMNANT SHOCKS 
Shock evolution 

A typical supernova explosion blows material into 
the interstellar/circumstellar medium (CSM/ISM) with 
velocities up to 10,000-20,000 km s^\ and a total kinetic 
energy of ~ 10^' erg. This ejected material is initial hot 
from the supernova explosion itself, but due to the high 
expansion rate, the ejecta cool very fast, even allowing 
dust to condense (as detected in the SNR Cassiopeia A 
in 18 1). The outer ejecta drive a shock wave through the 
CSM/ISM, heating the plasma to temperatures exceeding 
10^ K. 

The hot plasma resulting from the shock is at the out- 
side bound by unshocked CSM/ISM, but at the inside 
by fast moving, cool ejecta. The high pressure in the 
shocked medium starts driving a shock wave back into 
the ejecta, thereby heating the ejecta. This shock is called 
the reverse shock. Note that the direction of the reverse 
shock in the frame of the observer depends on the evo- 
lutionary phase of the supernova remnant: initially the 
reverse shock is also moving outward at high speed, and 
the faster moving ejecta "bump" into the reverse shock. 
At a later stage, occurring when the swept up CSM/ISM 
mass exceeds several times the ejecta mass, the reverse 
shock moves in an opposite direction to the forward 
shock, i.e. it moves toward the center. Once it reaches 
the center, all ejecta have been shock heated and the core 
of the SNR is hot. Such a SNR can no longer considered 
to be a young SNR, although some SNRs in this stage 
of the evolution still show evidence for highly metal en- 
riched ejecta [e.g. the LMC remnant Dem L71 .19..20,l . 

There are several analytic models for the shock struc- 
ture and evolution of SNRs |l2Tl|22l|23]. Fig.[T]shows the 
evolution of a SNR in a uniform ISM, as may be appro- 
priate for Type la SNRs, which are thought not to have 
a strong stellar wind prior to explosion. The figure il- 
lustrates the behavior of the reverse shock, whose radius 
peaks, in this particular case, at around 1000 yr 
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FIGURE 1. Shock radii and shock velocities as a function 
of time according to the Truelove-McKee n = 7 ,s = model 
|22|. The solid line represents the forward shock, whereas the 
dashed line represents the reverse shock. In the bottom panel 
the reverse shock velocity in the frame of the ejecta is shown as 
a dotted line. The dimensionless model was adjusted to fit the 
properties of Kepler's SNR |24|. 



The asymptotic behavior of the shock radius is that of 
the Sedov-Taylor self-similar evolution:i.e. : 

R=(2m—) , v. = — = --, (1) 
V po / at 5 t 

with E the explosion energy, and po the ISM density. 

The medium surrounding a SNR can be rather com- 
plex. In particular, around core collapse supernovae one 
expects several regions shaped by different evolutionary 
stages of the progenitor. On the main sequence a massive 
star blows a fast {v„ ^ 1000 kms^^), tenuous wind, but 
in the red supergiant phase a massive star loses mass with 
a rate of ^ 10""^ - 10"^ M©, in a slow (vv 10 kms"') 
dense wind. Finally, the most massive stars may become 
Wolf-Rayet stars, which again have a very fast, tenuous 
wind. As a result the CSM may consist of several layers 
separated by shells caused by the interactions of winds 
from different phases in the star's life. 

If one considers an SNR evolving in a red supergiant 
(RSG) phase, which may be quite common, one expects 
the following density profile ||2T1 : 

The young SNR Cas A may be currently in this phase 
Il25l |26I . The shock velocity in such a wind structure 
quickly follows the relation V, = |j, w hich is indeed 
close to what is found for Cas A ll27l[28l . 



Shock heating and temperature 
equiUbration 

The equations of mass-, momentum- and energy-flux 
conservation gives the foUowing relation between shock 
velocity, Vs, and temperature for a high Mach number 
shock [e.g.i29J: 

kT^^^mX- (3) 

T refers here to the average temperature, i.e. the av- 
erage kinetic energy per particle. However, the micro- 
physics of shock heating is not well known. SNR shocks 
are collisionless, which means that particle-particle col- 
lisions (Coulomb interactions) are rare, and insufficient 
to heat the plasma within a time comparable to the age of 
the SNR. Instead, the heating probably occurs through 
plasma waves. In such a case it is not quite known 
whether different particle species are heated to the same 
temperature. Alternatively, one may expect that all par- 
ticle species have the same velocity distribution. This 
could arise if plasma waves in the shock scatter the in- 
coming particles in different directions, with only small 
changes to the absolute velocities of the particles, which 
is Vs for all particles. In such a case the temperature 
for each particle species is different and scales with the 
mass, m,, of the particle: 

kTi=^m,V^. (4) 

In the post-shock region, full equilibration will take place 
on a time scale of log(net) ~ 12.5 lf30l . correspondingto 
~ 30,000 yr of a typical density of « 1 cm^^ (Fig.pb. 
Note that equilibration by particle-particle interactions 
depends on the charge of the particles and the mass ratio. 
For that reason electron-proton equilibration takes a long 
time, whereas proton-iron equilibration proceeds on a 
time scale of log(net) ~ 11-5. The parameter n^t is also 
important for the ionization process, and can be directly 
measured using X-ray line ratios}^ 

Electron temperature can be determined using X-ray 
spectroscopy, even employing CCD detectors. Ion tem- 
peratures can only be determined by measuring the ther- 
mal Doppler broadening of spectral lines. In order to 
eliminate Doppler broadening caused by line of sight 
motions due the SNR expansion, spectra should be ob- 
tained of the rim of a SNR, where only motions in the 
plane of the sky are to be expected. Most measurements 
of ion temperatures concern the hydrogen lines. Since 
hydrogen is quickly ionized, the hydrogen line emission 



arises from very close to the shock front. The presence 
of hydrogen lines is only possible, if some neutral hy- 
drogen is present in the ISM/CSM. The direct excitation 
in the post-shock gas gives rise to a narrow line, as the 
neutral hydrogen has not yet interacted with the shock 
heated plasma, and still has the temperature of the CSM. 
Charge exchange between an incoming neutral hydrogen 
atom and a proton in the heated gas, gives rise to a broad- 
ened line. The width of this broad component is a direct 
measurement of the proton temperature behind the shock 
[35 36). In the case of equilibration of temperatures, the 
proton temperature is equal to the mean plasma temper- 
ature. 

X-ray spectroscopy, in general, prodives a direct mea- 
sure of the electron temperature. The electron tempera- 
ture, together with rigt, determines the relative line ratios 
of lines of a given ion or atom. Also the characteristic 
cut-off in the bremsstrahlung continuum is directly re- 
lated to the electron temperature. In young SNRs, how- 
ever, synchrotron radiation is another source of contin- 
uum radiation, and the two mechanism are not always 
easy to tell apart. 

In principle, in X-rays one can also measure ther- 
mal Doppler broadening, but in practice the CCD de- 
tectors on board current X-ray observatories like Chan- 
dra and XMM-Newton pair great imaging capabilities 
with too poor spectral resolution to measure thermal 
Doppler widths. The high resolution grating spectrome- 
ters on board Chandra and XMM-Newton have a good 
spectral resolution, but they are slitless and the spec- 
tral quality is degraded for extended sources. Neverthe- 
less, X-ray Doppler broadening has been measured for 
SN1006 using XMM-Newton's Reflective Grating Spec- 
trometer li37J . This remnant has one of the lowest val- 
ues of rigt: rigtfv 2 x 10^ cm^^s. Optical, UV and X-ray 
spectroscopy of the northwestern part of this remnant 
reveals that ions, protons and electrons are out of tem- 
perature equilibrium ll38l[39l[37l . with the electrons be- 
ing much cooler than the oxygen ions (kTs w 1.5 keV, 
kTovu ~ 500 keV |f37l). 

Apart from SN1006, also in several other SNRs have 
ion temperatures been measured, but only using optical 
and UV lines |40, 41, 42, 43 1. Combining these mea- 
surements with electron temperatures obtained from X- 
ray spectroscopy shows that for the fastest shocks (Vs > 
1000 kms^') the electron temperature can be as low as 
10% of the proton temperature. For the slow shocks, like 
in the Cygnus Loop, electrons and ions seem to be in full 
equilibrium, with a turnover from equilibrium to non- 
equilibrium occuring around 600 km s^ ' |42 1 . Toward the 
end of this paper I will return to the issue of ion temper- 
atures in the context of cosmic ray acceleration. 



I will skip the details here, but for those interested, more information 
and background can be found in I30ll31ll32ll33ll34l . 
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FIGURE 2. Theoretical model of the collisional equilibra- 
tion process for a plane parallel shock model. It is assumed 
that the shock velocity is Vj = 4000 km s^ ' . The equilibration 
scales with the product of density and time, parameterized here 
as We'- The different lines indicate the temperature of different 
particle species. From bottom to top (on the left): electrons, 
protons, helium, oxygen, silicon, and iron. 



MULTIWAVELENGTH IMPRINTS OF 
EFFICIENT COSMIC RAY 
ACCELERATION 

The coming of age of Te V and X-ray astronomy over the 
last 15 years, has greatly increased our knowledge of cos- 
mic ray acceleration by SNRs. Here I will explain some 
of the observations that form the basis of our current un- 
derstanding. Not included here is TeV astronomy, which 
is amply addressed in other parts of these proceedings. 

Measuring magnetic fields using X-ray 
synchrotron rims 

The discovery by the ASCA satellite that the X-ray 
continuum emission from SN1006 is dominated by syn- 
chrotron emission Qj^has been the start of a number of 
discoveries in X-rays, which all relate to the cosmic ray 
acceleration properties of SNRs. 

One of these discoveries is that the X-ray synchrotron 
emission from SNRs is confined to a region very close 
to the shock front. In the SNRs Cas A, Kepler (SN1604) 
and Tycho (SN1572), the synchrotron emission comes 
from a region only a few arcseconds near the shock front 
BtI l48l . In fact, the identification of the X-ray syn- 
chrotron emission required the superior angular resolu- 



- This was anticipated by Reynolds & Chevalier 1461 . 



tion of Chandra. In some other remnants, SN1006, RCW 
86, RX J1713.7-3946, and RX J0852.0-4622, the X-ray 
synchrotron emission comes from a larger region, al- 
though some fine scale structure is present. 

The size of the X-ray emitting region is now generally 
believed to be determined by the magnetic field strength. 
The lifetime of a relativistic electron in a magnetic field 
is in cgs units: 

^ 637 

The typical photon energy Ep^ corresponding to an elec- 
tron energy E is: 

Ei,i, = 1.4E^Bj_keW. (6) 

The electrons are accelerated at the shock front by dif- 
fusive shock acceleration. At small scales their trajecto- 
ries are determined by scattering on plasma waves (dif- 
fusion), but on large scales their average motion follows 
the plasma. If the shock compression ratio is X ^ Po/Pu 
then the plasma velocity with respect to the shock is 
Vi — Vs/X- As the electrons are moving slowly away 
from the shock front, they lose energy, until they are 
no longer emitting X-ray radiation. Thus, the loss time 
corresponds to the width of the X-ray emitting region 
through: 

hoss = '^loss-Vs- (7) 

X 

By combining the measured width and the mea- 
sured/infeiTed shock velocity with the observed photon 
energy one arrives at an estimate of the post-shock 
magnetic field strength |7|, B, as illustrated in Fig. [3] 

Several groups |49, 8, 50 1 do not use the loss length 
scale, //oii, but instead assume that the length scale seen 
corresponds to the diffusion length scale, Idiff, i.e. the 
length scale at which diffusive motions are more impor- 
tant than the bulk plasma motion. However, in practice 
this method should give similar results, as long the syn- 
chrotron emission comes from electrons near the maxi- 
mum acceleration energy, where diffusive acceleration is 
balanced by radiative losses, since linff < lio^s for effi- 
cient acceleration ll34l ISTl l52l 

Magnetic field measurements of several SNRs have 
now been published, and indicate post-shock magnetic 
fields of 20-600 jxG |l2l[8l|53]|Sll9lll0l|55l|45l. This 
is higher than expected based on the compression of the 
average magnetic field in the ISM, B ^5 jJ-G, and, there- 
fore, suggests that some form of magnetic field amplifi- 



^ Interestingly, the two method have different assumptions: to estimate 
Idiff one has to assume Bohm-diffusion, whereas for estimating 
one has to assume a compression ratio. The fact that both methods 
give similar magnetic field values is an indication that both assumptions 
must be approximately valid 1341 . 
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FIGURE 3. A deep Chandra image of Cas A 1441 in the 4-6 keV continuum band (left). Note the thin filaments, marking the 
border of the remnant The remnant has a radius of about 2.5'. Right: The maximum cosmic ray electron energy versus magnetic 
field strength for the region just downstream of Cas A's shock front, as determined from the thickness of the filaments. The shaded 
area is excluded, because the filament width cannot be smaller than the minimum possible diffusion length [c.f.l?). (These figures 
were published before in \45i .) 
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FIGURE 4. Dependence of the post-shock magnetic field 
energy density S-/(8;rpo) on the shock velocity Vs- The dashed 
line shows a dependency and the dotted line a V'^ de- 
pendency. The input values can be found in |45|, except for 
SN1993J 1 56 1 and Kepler's SNR, for which the shock velocity 
in the southwest was taken from 1241 and the density from 1571 . 



cation mechanism is operating. Tiie most likely mecha- 
nism is non-linear growth of plasma waves induced by 
streaming of cosmic rays ll58l [TTl |59l . but some alter- 
native, perhaps complementary, mechanisms have been 
proposed jgOllgTl . 



Recently, high magnetic fields for SNRs were also 
reported based on small X-ray synchrotron brightness 
fluctuations in Cas A and RX J 171 3.7-3946 |62,„63,..64J. 
These findings support the presence of relatively high 
magnetic fields in young SNRs. 

It is not yet clear what the relation is between the 
ram pressure at the shock front, poV/, and the post- 
shock magnetic field pressure. Based on observations it 
has been argued that °^ PoV^ ||54l , whereas in ifTTI 
it has been argued that °^ Po^^- In the latter case, 
a proportionally larger fraction of the incoming kinetic 
energy is transferred to cosmic rays and magnetic fields 
for during the earliest life of an SNR, since then the 
shock velocity is highest. 

The problem with the current measurements of SNRs 
is that the dynamic range in density is quite high, but 
the dynamic range in shock velocity is small, since the 
magnetic fields can only be measured for X-ray syn- 
chrotron emitting remnants, and only SNRs with Vs > 
2000 km s^ ' |65 1 are expected to emit X-ray synchrotron 
radiation, whereas most Galactic remnants known have 
Vs < 5500 kms"'. However, for supernova SN1993J, 
which has a shock velocity ~ 20,000 kms^\ a high 
magnetic field of 64 G has been reported |56|. If one 
assumes that for this object the same mechanism is at 
work, one may try to distinguish between the two mag- 
netic field strength scalings. As Fig. [4] shows, includ- 
ing SN1993J favors a magnetic field strength scaling as 
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FIGURE 5. X-ray continuum spectral index map of Cas A 
based on tiie 1 Ms exposure by Chandra |66|. The gray scale 
coding indicates spectral number index with respect to F = 
—3.1. A hard spectral index (lighter pixels) is an indication for 
X-ray synchrotron emission. One sees here that hard spectral 
indices at the rims and in the central region, somewhat shifted 
to the west of the SNR. (Figure made by Eveline Helder, a color 
version can be found in L66J .) 

cx: P()V^, but one should regard this with some caution, 
since the magnetic field in SN1993J may have different 
origin than the magnetic fields in young SNRs. 

Magnetic field amplification at the reverse 
shock 

The high magnetic field of SN1993J raises the ques- 
tion how universal magnetic field amplification is. Does 
it require certain preconditions, like a medium strength 
magnetic field, to start with? 

Interestingly, it was recently reported that for the 
young SNR Cas A a large part of the X-ray synchrotron 
emission is coming from the reverse shock, in particular 
in the western region of Cas A 16611631 (Fig.|5]l. 

In the past, reverse shocks as sites of particle acceler- 
ation were often ignored (but see [67, ,68]), since it was 
assumed that the magnetic field in the ejecta was low, due 
to their large scale expansion. However, particle acceler- 
ation at the reverse shock, sufficiently fast to give rise to 
X-ray synchrotron emission, may be important for two 
reasons: 1) it shows that magnetic field amplification is 
rather universal, and that it does not require relatively 
large initial magnetic fields 2) some cosmic rays may be 
accelerated from metal rich material, which could lead to 
signatures in the cosmic ray composition. Indeed there 
is some evidence that 20% of the Galactic cosmic rays 
come from massive stars 1691 iTOl . 

As remarked by Luke Drury at the symposium, one 
should not overestimate the significance of the reverse 



shock for their contribution to the observed cosmic ray 
spectrum: cosmic ray acceleration is efficient when the 
physical shock velocity is high. For the reverse shock 
this occurs when the shock starts to move back toward 
the center of the SNR, until it reaches the center This is 
a relatively short period in the life of a SNR. Moreover, 
during this time the area spanned by the reverse shock 
is smaller than that of the forward shock, so the total 
number of particles entering the reverse shock is much 
less than the number of particles swept up by the forward 
shock. Both these aspects are illustrated in Fig. [T] For 
Cas A it is believed that the reverse shock velocity in 
the frame of the freely expanding ejecta is less than 
2000 kms^'|23|, except in the western region, where 
somehow the reverse shock is almost at a standstill in 
our frame, and the ejecta are shocked with the ejecta free 
expansion velocity of 4000 - 8000 km s~^l66l. 

What is the maximum energy cosmic rays 
can be accelerated to? 

Several authors have pointed out that the high mag- 
netic fields that have been inferred from X-ray syn- 
chrotron emitting rims indeed allow protons to be accel- 
erated to > lO'^ eV, and heavy ions can in principle ac- 
celerated to even higher energies [e.g. 34, 51 1. It is also 
clear that the magnetic field energy density scales with: 

B^-poV", (8) 

with a = 2 or 3. As a consequence higher magnetic fields 
are present early in the life of a SNR when the shock 
velocity is higher This has led to the suggestion that the 
highest energy cosmic rays are accelerated early on, and, 
as the magnetic field drops, those particles escape first. 
The peak energy is thus a function of time. The observed 
cosmic ray spectrum is in such a scenario a superposition 
of cosmic rays released over an extended period of time 

m 

Note that not only does the magnetic field energy de- 
pend on Vs, but also on po. This means that the highest 
energies are reached for shocks in dense environments. 
Some supernovae have by their very nature high circum- 
stellar densities. As indicated by Eq.[2] the density around 
a massive, wind blowing star, is highest for slow wind 
speeds, v,,,, and high mass loss rates. So it is likely that 
the highest cosmic ray energies are obtained by SNRs 
developing in a RSG wind. An example, as mentioned 
above, is Cas A. In fact, an optical spectrum of the light 
echo of the supernova explosion reveals that it is the rem- 
nant of a Type lib supernova, very similar to the bright 
radio supernova SN1993J [72|! Cas A's progenitor prob- 
ably had only a very short Wolf-Rayet star phase, if any 
atall||25j|73l|26l. 



An SNR like Cas A reaches the self-similar shock 
evolution very early on, in which case °= f^^^ and 
-1/3 lEl- One can use this similarity evolution 
to estimate the maximum proton cosmic ray energy as 
a function of time: First note that BiV^, with 

tacc the acceleration time, for which we can take the age 
of the SNR. Using Eq.|2]and Eq.[8]one finds: 

Bl oc J-f-^/^ocf-^/^""/*'. (9) 



We therefore find for the maximum proton energy in a 
young SNR evolving in a dense stellar wind: 

E oc ,i-4/3-a/6 



for a = 3 this means E,,,,^,^ °^ t ^1^, 

^-2/3 



(10) 

and for a = 2 this is 

Emax °^ f 

The value for the magnetic field found for Cas A (Z? w 
0.5 mG) is sufficient for E^ax presently to be ~ lO'^ eV. 
However, Cas A's ability to accelerate cosmic rays was 
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probably even better in the past: according to Eq 
and using a = 3, when Cas A was only 30 yr old the 
maximum proton energy was £■„,£„ 7 x lO'^ eV and its 
post-shock magnetic field strength 7 mG. Note that the 
evolution of the magnetic field ensures that the diffusion 
length scale (assuming Bohm diffusion) does not exceed 
the size of the remnant, since B evolves faster with t 
than Rs- The density profile makes that the flux of 
particles does not depend on the shock radius. 

All this implies that probably most of the Galactic cos- 
mic rays with E > 10'^ eV were accelerated by Type 
II/Type lib supernova remnants, during the first 100 yr 
of their lives. Most massive stars explode while in their 
RSG phase. Only stars more massive than 25 M0 (a mi- 
nority), probably explode in the Wolf-Rayet stars phase. 
Their SNRs are probably less capable of acceleration 
cosmic rays to very high energies, unless they explode 
with larger energies, as in the case of hypernovae. 



Measuring compression ratios using X-ray 
imaging 

Magnetic field amplification is probably a result of 
cosmic ray streaming. The presence of relatively large 
magnetic fields in young SNRs, therefore, implies by it- 
self efficient cosmic ray acceleration. I use the word "ef- 
ficient" here in contrast to fast acceleration. Fast accel- 
eration means that cosmic rays are accelerated to high 
energies, efficient that the cosmic ray content of SNRs is 
energetically important. 

Efficient cosmic ray acceleration affects both the evo- 
lution of the SNR itself ll74l . and the cosmic ray spec- 
trum 1751 1761 1771. If cosmic ray acceleration is very ef- 
ficient the cosmic rays provide back-pressure to the un- 
shocked plasma, which gives rise to a concave spectrum. 



i.e. it is steeper than the test particle spectrum at low 
energies and flatter at high energies. If the internal en- 
ergy of the post-shock gas is dominated by relativistic 
cosmic rays the compression ratio increases. A standard 
high Mach number shock in a monatomic gas (7= 5/3) 
will have a compression ratio of % = (7+ 1)/(7^ 1) = 4, 
whereas for a gas dominated by relativistic particles (7 = 
4/3) this IS, X — 1 ■ Compression ratios in excess of 7 are 
possible, if energy losses are taken into account. 

In old SNRs radiative energy losses lead to 
higher compression ratio, since shock velocities of 
Vs < 200 kms^^ produce plasma with a temperature 
of < 10^ keV, near the peak of the cooling curve. This 
gives rise to the formation of the filamentary structures 
radiating in forbidden lines, which make for beautiful 
Hubble Space Telescope images, of SNRs like Vela, the 
Cygnus Loop or, N49. 

For young SNRs, which efficiently accelerate cosmic 
rays, energy may leak out of the plasma due to the escape 
of cosmic rays. The highest energy cosmic rays are the 
most likely ones to escape. Escaping cosmic rays are 
only energetically important for flat, or concave spectra. 
In that case the highest energy cosmic rays may contain 
a significant fraction of the total internal plasma energy. 
Although concave spectra are expected in most efficient 
cosmic ray acceleration models, this is by no means a 
certainty yet; see the contribution by V. Zirakashvili in 
these proceedings. 

Higher compression ratios due to cosmic ray escape 
will give rise to different ratios of the forward and re- 
verse shock radii or of the forward shock and contact dis- 
continuity radii |74|, the contact discontinuity being the 
boundary between shocked supernova ejecta and shock 
heated CSM. 

Indeed, for both Tycho's SNR |55 1 and SN1006 [TS] it 
has been reported that the contact discontinuity lies very 
close to the forward shock. SN1006, a newly detected 
H.E.S.S. source, is an interesting SNR in this respect, 
since the X-ray synchrotron emission, and also most of 
the radio emission, is confined to the northeastern and 
southwestern regions. If we take the shocks in these 
regions sites of efficient cosmic ray acceleration, one 
may expect that in those regions the contact discontinuity 
and forward shock are closer to each other than in the rest 
of the remnant, due to a higher compression ratio. This 
is indeed what is found, but surprisingly also in the other 
regions of SN1006 the contact discontinuity is closer to 
the forward shock than expected |78 1. 

This suggests that apart from a higher compression 
ratio also some other mechanisms (e.g. ejecta clumping, 
hydrodynamical instabilities) are at play in bringing the 
ejecta close to the shock front. So SNR morphology hints 
at high compression ratios, signifying efficient cosmic 
ray acceleration, but a quantitative results cannot yet be 
obtained. 



Radio polarimetry and magnetic field 
ampUfication 



High resolution spectroscopy and cosmic 
ray acceleration efficiency 



Fast cosmic ray acceleration requires a small diffusion 
constant, for which, apart from a relatively high mag- 
netic field, also a turbulent magnetic field is required, 
i.e. 5B/B ^ 1. Magnetic field amplification by growth 
of plasma waves due to streaming instabilities ifTTl . natu- 
rally provides both the magnetic field strength and its tur- 
bulence. Observationally, turbulent magnetic fields are 
implied by the fact that near the maximum electron ener- 
gies the diffusion length scale and advection (loss time) 
length scales are similar (see p4| and above). 

Other evidence that the magnetic fields are turbulent 
comes from the long known observations of radio po- 
larization of SNRs. These indicate that old SNRs have 
preferentially tangential magnetic fields, whereas young 
SNRs have radial magnetic fields [79J. Some recent re- 
sults regarding magnetic field structure can be found for 
Cas A pTl, Kepler's SNR |80|, RCW 86 fST] and the old 
SNRPKS 1209-51 |82|. 

The tangential magnetic fields in old SNRs are caused 
by the compression of the ISM magnetic fields, which 
have some large scale coherence (see the contribution 
by R. Beck in this volume). Shocks only compresses the 
tangential component of the field. In those regions where 
the uncompressed magnetic fields are perpendicular to 
the shock normal, the post-shock magnetic field will 
be enhanced, increasing the radio brigthness of these 
regions. 

The radial magnetic fields in young SNRs are not 
that easily explained by MHD simulations [83] l84l . 
Rayleigh-Taylor instabilities at the contact discontinu- 
ity, indeed, produce radial magnetic fields f83). How- 
ever, near the shock front, simulations produce tangen- 
tial, not radial magnetic fields. Efficient cosmic ray ac- 
celeration, may indeed be the missing ingredient. This 
can happen in two ways. First of all, efficient cosmic 
ray acceleration gives rise to turbulent magnetic fields. 
Since these have no preferred direction, strong tangen- 
tial fields where the shock normal is perpendicular to 
the magnetic fields are avoided |84|. Secondly, enhanced 
shock compression will bring the Rayleigh-Taylor insta- 
bilities closer to the shock front ||85]| . However, it is not 
clear whether the compression is sufficiently enhanced 
to indeed produce radial magnetic fields close to the for- 
ward shock. Zirakashvili reported at this symposium that 
MHD simulations incorporating magnetic field ampli- 
fication tends to make the downstream magnetic field 
stretch in a radial direction ll59l . 

It, therefore, seems likely that the observed radial 
magnetic field structure in young SNRs is either facili- 
tated or a direct by-product of efficient cosmic ray accel- 
eration and magnetic field amplification. 



In the literature on thermal Doppler broadening of 
lines from SNRs, one has concentrated mostly on the 
equilibration of electrons and ions, and the consequences 
for deriving shock speeds using Eq. |3] or Eq. |4] Rarely 
has the discussion focussed on the influence of cosmic 
rays on the measured temperatures. However, cosmic 
rays may contain a large part of the internal energy, and 
change the adiabatic index of the plasma. These effects 
are not taken into account in Eqs.[3] |4] which are usually 
used to derive shock parameters. 

An exception is the case of the Small Magellanic 
Cloud remnant IE 0102.2-7219. As reported in [86J the 
electron temperature of this remnant is too cold com- 
pared to the measured shock velocity, even allowing for 
non-equilibration of electron and ion temperatures. A di- 
rect measurement of the ion temperature in this remnant 
is difficult, as this oxygen rich remnant does not emit op- 
tical hydrogen lines. 

How much does efficient cosmic ray acceleration in- 
fluence the post-shock plasma temperatures? The answer 
depends on two parameters 1) the fraction of the pressure 
contributed by non-thermal components, 2) the energy 
escape fraction associated with cosmic rays. For a strong 
shock one can parameterize these with: 



Pnt 



Pt+Pnt 



(11) 



the relative pressure of the non-thermal components 
(cosmic rays and magnetic fields) [c.f. "871 iTTll . and the 
shock compression ratios, which for a strong shock is: 



G+VG2-(1-£,,,)(2G-1) 
^ 1-e... ' ^ ^ 

with Eesc the fraction of the incoming energy flux that 
is taken away from the shock by cosmic rays, and G = 
5W+ ll^One can apply these expressions to the follow- 
ing relation between plasma temperature and shock tem- 
perature: 



kT. 



(l-w)-(l 



(13) 



which is a generalized version of Eq.Hl as can be seen by 



inserting w = and X=^- El- 13 andfOare graphically 



depicted as a function of w and e^^c in Fig |6] 

It is clear from these equations that high resolution 
spectroscopy offers the opportunity to learn about the in- 
ternal energy budget of the plasma. Measuring kT-p and Vs 



^ I assume here that the non-thermal contributions to the pressure have 
an adiabatic index of 7 = 4/3. 
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FIGURE 6. Graphical representation of the effects of cosmic ray escape and relative cosmic ray pressure (w) on the compression 
ratio (left) and the post-shock temperature (right). In the right hand figure, the values are given as correction factors with respect to 
the standard strong shock expressions (Eq. |3|4^ . 
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FIGURE 7. Hce line profile of the northeastern part of RCW 
86, as observed with the VLT (Helder et al. in preparation). 
The line consists of a narrow peak, from direct excitation, and 
broad wings caused by charge exchange between the incoming 
neutral hydrogen atoms and the shock heated protons. The 
FWHM of the broad wing is a direct measure of the proton 
temperature. The narrow peak is broadend due to the relatively 
broad line spread function (LSF) of 330 kms^'. Preliminary 
results give a FWHM of (1270 ±60) kms"', after correction 
for the LSF. (Figure provided by E. Helder) 



So far the interpretation of most measurements of kT 
using optical, UV and X-ray spectroscopy have ignored 
the role of cosmic ray physics (but see |88|). To some 
extent this seemed unnecessary, because spectra were 
taken from locations from which no strong X-ray syn- 
chrotron radiation is emitted, like the northwestern re- 
gion of SN1006. The northwestern region of SN1006 has 
bright Ha emission, which made it easy to obtain high 
quality spectra. The X-ray synchrotron regions, on the 
other hand, only show weak Ha emission. That is unfor- 
tunate, since at these synchrotron rims cosmic ray accel- 
eration is likely to be more efficient. It is not quite clear 
whether the lack of Ha emission from the X-ray syn- 
chrotron rims is a coincidence. Three possible reasons 
for this anti-correlations are 1) the presence of neutral 
hydrogen damps plasma wave, decreasing the efficiency 
for cosmic ray acceleration |89|; 2) cosmic rays diffus- 
ing away from the shock ionize the CSM, supressing the 
number of neutral hydrogen atoms entering the shock; 3) 
bright Ha emission comes from the densest regions, but 
as a result of the higher density the shock velocity has de- 
celerated, thereby diminishing the acceleration efficiency 
f52l. 



Knot g in Tycho 's SNR 



independently gives a handle on cosmic ray acceleration 
efficiencies. This should be combined with a measure- 
ment of kTfi in order to test for ion-electron equilibration. 
Ideally, one would like to measure also the compression 
factor independently, in order to alleviate the degeneracy 
between w and Eesc- On the other hand, it seems unlikely 
that w < Eesc, which constrains part of the degeneracy. 



A case in point is the Ha-bright "Knot g" in Tycho's 
SNR. The broad Ha line component gives a proton tem- 
perature of 5.9 ± 0.9 keV, indicating a shock velocity of 
1700-2200 kms^', depending on whether ion-electron 
equilibration is assumed or not. The average shock veloc- 
ity of Tycho's SNR is 3 100 km s^ ' , based on radio proper 
motion studies li90l . and assuming a distance of 2.5 kpc. 



Without cosmic ray acceleration the expected tempera- 
ture is kTp = 19 keV for an unequilibrated plasma, or 
kT — II keV for an equilibrated plasma. However, for 
"knot g" the shock velocity is probably much lower, 
since the radio proper motion at this location indicates 
an expansion parameter Vj (Rt/t) ~ 0.25, rather than the 
0.46 found in the rest of the remnant. The reason is prob- 
ably that the blast wave at "knot g" is encountering dense 
material slowing down the shock, but also resulting in 
bright Ha emission. For the lower expansion velocity of 
"knot g" the measured plasma temperature is consistent 
with a shock without cosmic ray modification. This is at 
odds with the above mentioned findings that all over this 
remnant the contact discontinuity is too close to the for- 
ward shock as a result of cosmic ray acceleration ||55]| : 
cosmic rays do not seem to contribute much to the inter- 
nal energy in "knot g". 



The low plasma temperature in LMC SNR 0509-67.5 

One of the fastest expanding young SNRs is 0509- 
67.5 in the Large Magellanic Cloud. This remnant has 
an age of about 400 yr ll45l |9T1 and high resolution X- 
ray spectroscopy with XMM-Newton's RGS experiment 
shows line broadenings of (7,. w 5000 kms^^, which 
is dominated by kinematic Doppler broadening, as the 
whole SNR is observed |92|. The broadening provides 
a lower limit to the shock velocity, since the plasma ve- 
locity directly behind the shock is 3/4^5, and most of 
the line emission comes from the reverse shock region, 
which moves probably 30% slower than the plasma im- 
mediate downstream of the shock velocity. This means 
that the shock velocity is at least « 6700 kms^\ but 
likely larger than Vs ~ 8000 kms^'. The hydrogen line 
broadening has been measured to be 3710±400 kms^^ 
(FWHM) 1,43 J , implying a proton temperature of kTp = 
26 keV. This translates into a shock velocity 3000- 
5400 km s^^ which is much lower than the actual shock 
velocity. The ratio of the measured temperature and the 
expected temperature is ^ 0.2 (for Vs « 8000 kms^'), 
indicating, according to Fig.|6] either a high partial pres- 
sure not coming from the thermal plasma, w — 0.7, or 
a combination of lower w with a relatively high cosmic 
ray escape fraction. Our knowledge of cosmic ray accel- 
eration in 0509-67.5 would be more complete, if either 
X-ray synchrotron emission or TeV 7-rays would be de- 
tected. Unfortunately, for the moment the remnant is too 
distant to identify X-ray synchrotron emission from the 
shock front. 



The plasma temperature in an X-ray synchrotron 
emitting region in RCW 86 

Ideally one would like to measure post shock temper- 
atures of region that emit both X-ray synchrotron radia- 
tion and/or TeV 7-rays. In the future this can be under- 
taken with the next generation of high resolution, imag- 
ing X-ray spectrometers, as are planned for the Japanese 
NEXT mission or the ESA/NASA/JAXA International 
X-ray Observatory (IXO). X-ray and UV spectroscopy 
have advantages over Ha spectroscopy, since X-ray line 
emission does not require a partially neutral pre-shock 
gas. 

In the mean time one has to concentrate on the faint 
Ha line emission coinciding with X-ray synchrotron 
emission. In a new study Helder et al. (in preperation) 
have done that for the SNR RCW 86, the possible rem- 
nant of AD 185 II93. RCW 86 has recently been detected 
in TeV 7-rays by H.E.S.S. (Hoppe et al. these proceed- 
ings). Its X-ray emission is a mix of X-ray synchrotron 
, and thermal radiation. The X-ray emission indicates 
strong density contrasts, and probably also strong gra- 
dients in shock velocities along the shell L52J . Ha emis- 
sion has been detected all around the shell ||94l , but the 
Ha emission is weak in the northeastern side, where, the 
emission is dominated by X-ray synchrotron radiation. 
As a consequence ESO's Very Large Telescope is needed 
to obtain Ha spectroscopy of this region (Fig. [7] shows 
the prelimanary spectrum). Surprisingly, the line width 
is 1270 km s^'(FWHM), corresponding to a proton tem- 
perature of kTp — 3 keV, (ignoring proton-electron equili- 
brartion). This is much larger than in other parts of RCW 
86, where typically line width of 500 kms^We found 
|41|, but it is smaller than expected for an X-ray syn- 
chrotron emitting shock for which Vs > 2000 km s^ ' ll65l . 
The plasma temperature appears to be too low, which 
can be explained, if cosmic rays have taken up part 
of the shock energy. For a reasonable value of Vs w 
2500 kms^'. Fig. |6j suggests w = 0.6 for an escape frac- 
tion of Eesc — 0. 

However, an alternative scenario to explain the low 
plasma temperature in RCW 86 is that the X-ray syn- 
chrotron emitting electrons were accelerated in the past, 
when the shock velocity was higher. This requires a long 
loss time for the electrons, corresponding to a relatively 
low magnetic field. Indeed, the value for the magnetic 
field derived from the X-ray emitting region is 24/iG 
|52|, corresponding to a sufficiently long loss time of 
^ 500 yr The new H.E.S.S. results are in agreement with 
such a low magnetic field, provided that the TeV emis- 
sion has a leptonic, rather than a hadronic origin. 

So, unfortunately, no clear cut answer can be given 
about the cosmic ray content of RCW 86, as long as 
the shock velocity is not directly measured, or the 7-ray 



emission process is identified. However, RCW 86 illus- 
trates how seemingly remote areas of astrophysics, 7-ray 
astronomy and optical spectroscopy, are both needed to 
obtain answers about cosmic ray acceleration efficien- 
cies. 



SUMMARY AND CONCLUSION 

The topic of these proceedings is 7-ray astrophysics, a 
field of astrophysics that is still relatively young, but has 
already changed changed our perception of cosmic ray 
acceleration in SNRs. However, as I have indicated here, 
more traditional fields of astronomy valuable clues about 
cosmic ray acceleration in young SNRs, as well. 

In reality, not all signatures of cosmic ray acceleration 
are unambiguous. For example, the evidence for high 
compression ratios is tantalizing |55 78 1, but the high 
compression ratios all over the remnants, even those 
region where cosmic ray acceleration seems not to be 
efficient, is confusing, and requires further investigation. 
However, it is important that we know now for what kind 
of multiwavelength signatures of cosmic ray acceleration 
to look, and what kind of observations are needed to 
make further progress. 

Comparing our current knowledge of cosmic ray ac- 
celeration by young SNRs with what was known a 
decade ago, one sees that a lot of observational progress 
has been made. My guess is that most cosmic ray physi- 
cists are a much more optimistic about the idea that 
young SNRs are the sources of cosmic rays up to the 
"knee" than 15 years ago. 
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